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1 Introduction 



In this talk, I would like to speak of the problem of the stability of three- and four-body non- 
relativistic purely Coulombic systems. A system will be said to be stable if its energy is lower 
than the energy of any subdivision in subsystems. This is a restrictive definition of stability, 
because besides that there are other useful notions: "metastability" and "quasi-stability" on 
which we shall say only a few words later. 

The works I will present are due, in what concerns the three-body case, to J.-M. Richard, 
T.T. Wu and myself. The four-body work is due to J.-M. Richard in collaboration with various 
other persons (including J. Frohlich !) 

The reasons why I decided to choose this subject for this workshop dedicated to Walter 
Thirring are that I know that Walter is interested in that topic and that the tools which are 
used are found precisely in Walter's celebrated quantum mechanics course: concavity, scaling, 
and the Feynman-Hellmann theorem. One tries to avoid, as much as possible, numerical cal- 
culations, or to use already existing numerical calculations in particular cases. I will speak 
of: 

i) three-body systems with equal absolute value of the charge, i.e., —e + e + e, or +e — e — e, 
since it is clear that +e + e + e is unbound. Then binding or no binding will depend on 
the masses; 

ii) three-body systems with unequal charges; 

iii) four-body systems with charges with equal absolute value. It will be mostly +e + e — e — e. 
However, I shall say a word on +e — e — e — e. 

2 Three-body case: equal | charges | 

The problem we discuss now is whether a system of three charged particles (1,2,3), 1 having 
charge +e and 2 and 3 charges — e, is stable or will dissociate into a two-body system and an 
isolated particle, (l,2)+3 or (l,3)+2. The system will be stable if the algebraic binding energy 
of the (123) system is strictly less than the binding energy of both (1,2) and (1,3). If, on the 
other hand, the infinimum of the spectrum of the (123) system coincides with the lowest of the 
(12) and (13) binding energies the system will be unstable. 

This is an old problem which has been treated in many particular cases. For instance, long 
ago, Bethe has shown that the hydrogen negative ion (pe~e~) has one bound state and Hill 
has shown that there is only one such bound state with natural parity J2|, and Drake has also 
shown that there exists an unnatural parity state [[3] and finally Grosse and Pittner M have 
shown also that this unnatural parity state is unique. In what follows we shall treat only the 
natural parity states, i.e., states such that P = (— 1) L , where L is the total orbital angular 
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momentum (we neglect spin interactions!). For three particles there is no problem with the 
Pauli principle even if two of them are identical fermions, since we can adjust the spin. 

Wheeler || has also shown that the system e + e~e~ is bound, and, more generally, Hill 
has shown that any three-body system in which the two particles with the same sign of the 
charge have the same mass is stable. This covers the two previous cases. 

As an example of an unstable system (there are many others !), we can give the proton- 
electron-negative muon system, for which a heuristic proof was given by Wightman in his thesis 
|| and a rigorous proof was given by a collaboration including Walter Thirring himself [^]. 

Richard, Wu and myself |§ have tried to organise the results on stability, and, by using 
simple properties, save numerical calculations. From the reactions we had from experts on 
numerical calculations we believe that this was not totally useless. The three-body Schrodinger 
equation reads 



' - -^—A 2 ip - -^—A 3 ip + 



2mi 



2mo 



2m 3 



e 2 e 2 ^ e 2 



ri3 r 23 



(1) 



and the corresponding two-body equations can be obtained by omitting some terms. It is 
obvious that we have scaling properties: 

i) the charges can be multiplied by some arbitrary number without changing the stability 
problem; 

ii) the masses can also be multiplied by an arbitrary number, so that the stability problem 
depends only on the ratio of the masses, i.e., of 2 parameters. 



It will be convenient to introduce some variables: 



the inverse of the masses 

111 

xi = — x 2 = — x 3 = — (2) 
rrii 1^2 1^3 

then the ground state energy of the system will be concave in x%, x 2 , x 3 , and, in particular, 
concave in x\ when x 2 and x 3 are fixed (and circular permutations!); 

the constrained inverse of the masses 

cti = etc. (3) 

X\ + x 2 + x 3 

such that 

cti + a 2 + a 3 = 1 (4) 
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With these new variables, any system of three particles can be represented by a point in a 
triangle, ai,a 2 ,a 3 being the distances to the sides of the triangle. Figure 1 represents such a 
triangle with a few points representing some three-body systems. 




Figure 1: 



It is of course sufficient, for the time being (i.e., for equal charges of 2 and 3) to consider 
the left half of the triangle, i.e., to assume vn 2 > 777-3. Let us remember that since we have, 
according to Hill's theorem, strict stability for m 2 = m 3 , i.e., a 2 = «3, there will be some 
neighbourhood of the line a 2 = a 3 where we shall have stability. However, not all systems will 
be stable. We have already mentioned the pe~fi~ system as unstable. Another point where 
instability is obvious is the left summit marked 3, where we have two infinitely heavy particles 
with opposite charge producing zero attraction on the third particle. There is, therefore, an 
instability region in the left half triangle. 

We have proved three theorems on the instability region in the left half triangle. 

Theorem I 

The instability region in the left half triangle is star-shaped with respect to summit 3. 

The proof is based on the Feynman-Hellmann theorem combined with scaling, take a point 
P (Fig. 2) where the system is unstable or at the limit of stability. First we use the variables 
xi,x 2 ,x 3 . From the Feynman-Hellmann theorem, d ^ 23 - ) > 0, if x\ and x 2 are fixed. The 
binding energy of the subsystem 12 is fixed. Hence the residual binding can only increase 
(algebraically). £3 moves from x 3 (P) to infinity. The image of this in the rescaled a variable is 
the segment, P3, where 0:1/0:2 = constant. If there is no binding at P there is no binding on 
the whole segment. 

Theorem II 

In the left-half triangle, the instability region is convex . 

Take two points P' and P" on the border of the stability domain inside the triangle with the 
a variables. At P' and P" we have E P ,{12) = £ P /(123), E P »(12) = £ P »(i23). It is possible to 
find a linear rescaling P -> M such that E M ,{\2) = E M „(12) = E M ,(123) = E M »(123). Then 
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Figure 2: 



one can interpolate linearly between M' and M": 



M\ = AM' + (1 - \)M" , 



< A < 1 



For any Ma, £m a (12) = const, and £a/ a (123) is concave in A, and therefore E Mx (123) > 
£jv/'(123) = Em>'(123)- Returning to the original variables a\Ct<iOi?, and noticing that the scaling 
is linear we see that, on P'P" we have £(123) > £'(12). Hence we have instability (Fig. 2). 

There is, in fact, a more refined theorem, which we found, following a question by the late 
V.N. Gribov during a seminar in Budapest in 1996. 

Theorem III 

The domain (in the left-half triangle) where 



The meaning of this theorem is that the lines along which the relative binding is constant 
have a definite convexity. Note the sign of the inequality because £(123) and £(12) are both 
negative. We believe that the proof is essentially obvious, since, in the previous theorem, one 
goes through a rescaling, replacing P and P' by M and M' where the two-body energies are 
equal. Theorem II is of course becoming a special case of Theorem III, with e — > 0. The dotted 
line on Fig. 2 corresponds to some positive value of e. 

Let us give a very simple application of Theorem II. We know that, according to Glaser et 
al., the system p OQ A~e~ is unstable if m^- > 1.57m e - (p^ means a proton with infinite mass). 
Similarly, we know that, from the work of Armour and Schrader || the system p OQ A + B~ is 
unstable if m^/mg- < 1.51. 

This means that pooe~e + is unstable (not because of annihilation that we neglect, but of 
dissociation into Pooe~ and e + ). In Fig. 3, p OQ A + B^ and pooA~e~ with the limit masses 
correspond respectively to X and Y. Any point to the left of the segment XY corresponds, 



£(123) 
£(12) 



< 1 + e , 



e > 



is convex. 
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according to Theorem II to an unstable system. Therefore, pe + e , , with the actual mass 

of the proton, are unstable, and one can go up to pz~z + which will be unstable if m p /m z > 2.2. 




Figure 3: 



We obtain too that the system p/j, e , with the actual proton mass is unstable, and also 
(disregarding again annihilation) ppe~ , pp\i~ . 

One can also use convexity to get results in the opposite direction, i.e., prove that certain 
three-body systems are stable. We know that the systems represented by a point on the vertical 
bissector of the triangle are stable. In practice we know more than that, namely we have an 
estimate or more exactly a lower bound of the absolute value of the binding energy of many 
systems by using variational calculations and, in fact, by playing with convexity again it is 
possible to have a lower bound of the absolute value of the binding energy at any point on the 
bissector which corresponds to a 2 — a 3 , < a 2 < 1. 

Now we use convexity along a horizontal line a\ = const. 

The systems ai,a 2 ,a 3 , a±, a 3 ,a 2 represented in Fig. 2 by Q and Q', are of course com- 
pletely equivalent. Hence 

i^!23(«l, «2, «3j = -&123(ttl, «3, OL 2 ) < -&123 I «1, ^ , ^ I , 

by convexity and 

/ a 2 + a 3 a 2 + a 3 \ ( \ ( a 2 + a 3 \ ( \ / 1 - on 

Ei 23 lai, , I = I l + g(ai) )E 12 \ a 1: I = \ l + g(a 1 ) )E 12 \ a 1: — - — 

where g represents the relative excess in binding energy. We are assured of stability if 

E 12 (a!,a 2 )> (l + g(a 1 )^jE 12 [a 1 , 1 

i.e., if 

e 2 1 / . ,\e 2 2 

— < 1 + ^("i) — 

2 ai + a 2 V yV V 2 ! + 
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In this way, it is possible to prove that the system pdfi~ } important for fusion processes, is 
stable, through it is off the diagonal. One would like to show also in this way that 7i + fi~fi + and 
/i + 7T + 7r~ are stable, but these considerations are not sufficient. There is a hint that they are 
stable because, from explicit calculations at ot\ = one sees that this method tends to give a 
band of stability which is two times narrower than the real one and this is just what one needs. 



3 Three-body case. Unequal charges 

On this topic, Richard, Wu and myself have published one paper |l(| and one in preparation 



of which I shall give some of the really new results. We have the right to take q\ — 1, the charge 
of the particle which is opposite to the other two, of the same sign, q 2 and q^. 

A) Unequal charges, but q 2 = q% 

This is the simplest case, very similar to the case of all equal charges. For fixed q 2 — 93 we 
can again represent a system with the variables a\a 2 ot^, and, on the bissector of summit 1, the 
energies of the subsystems (12) and (13) are equal. The fact that the instability regions are 
star-shaped with respect to 3 for the left-half of the triangle and to 2 for the right-half persists, 
and as well the convexity of the instability regions. There are two major differences which are: 



i) that if q2 = qs < 1, all three-body systems are stable, because near summit 3, for instance, 
the subsystem (12) is very compact and exerts a Coulomb attraction at long distances 
on particle 3; it may seem strange that as q 2 — > 1 part of the triangle becomes unstable, 
but this is just due to the fact that the binding energy, in that region, tends to zero as 
<?2 -> 1; 

ii) that if q 2 = q% is large enough, stability disappears completely. 



Figure 4 summarizes the situation. For q 2 > 1 but very close to 1, there is no qualitative 
difference but for a certain critical value 1 < q 2c < 11, the stability band breaks into two 
pieces, and from calculations by Hill and collaborators [|11|] stability near a 2 = a 3 = 1/2, 
ax = disappears completely for q 2 > 1.1, and from the calculations of Hogreve it disappears 
near a 2 = a 3 = 0, ot\ = 1 for q 2 > 1.24. From convexity, it hence disappears completely along 
the segment joining a 2 = = and a 2 = = 1/2, and from the star-shaped property, there 
is no stability at any point in the triangle for q 2 > 1.24. 
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q2< 1 



P2=1 



$2 ~ 1 





q2 - 1 .09 



1.1 <q 2 <1.24 




q 2 > 1 .24 



Figure 4: 



B) Unequal charges, but q 2 ^ q% fixed 

First we continue to use the a±, a 2 , ot% variables to describe the three-body system for fixed 
charges. A fundamental difference is that the bissector of summit 1 of the triangle no longer 
plays a special role. It is, instead, the line along which Eyi = Ei 3 , i.e., 



it 
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a\ + a 2 on + a 3 



or q\ (l - a 2 ^j = q\ (\ - a 3 ) , 



which becomes important. This line goes through the point a 2 — (X3 — 1, symmetric of summit 
1 with respect to the line a± = (Fig. 5). The line divides the triangle into two subregions. If 
we decide to take q 2 > q 3 , E 12 < E± 3 in the left region which contains the summit 1. 




Figure 5: 



If q 2 and ^3 are both less than 1, we have again stability everywhere. If q 2 > 1 with 
^3 < 1, part of the triangle becomes unstable. Various scenarios are shown in Fig. 6. Notice 
that summit 3, on the left is unstable, together with some neighbourhood. This is because 
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q 3 < q 2 < 1 



q 3 < q 2 = 1 
q 3 not too small 




q 3 < q 2 = 1 q 3 < 1 < q 2 q 3 < 1 < q 2 

q 3 small q 3 /q 2 not too small q 3 /q 2 small 



Figure 6: 

Oil = a 2 = corresponds to a very compact system with either very weak rapidly decreasing 
attraction on particle 3 (if q 2 — 1) or repulsion (if g 2 > 1)- 

For q 2 and q 3 both > 1, we have only a partial understanding of the situation. If they are 
both not too large, a neighbourhood of the dividing line will survive as stability region plus a 
neighbourhood of summit 1. Summits 2 and 3 will definitely be unstable. 

We believe that if q2 and #3 are sufficiently large there is no stability at all, for any mass, 
but this is very difficult to implement quantitatively except in two places: 



1. at summit 1, which corresponds to a 2 = a 3 = 0, and is the Born-Oppenheimer limit, for 
which Hogreve [O has shown that one has instability if either q 2 or q 3 > 1.24; 



2. at the point «i = 0, a 2 
there is no stability. 



a 3 



1/2 where Lieb's theorem applies: if l/q 2 + l/g 3 < 1 



This implies by convexity that if q 2 and q 3 > 2 one has instability for a 2 = 03, and by the 
star property for a 2 < 0:3 (if <?2 > ffe)- For a considerable improvement, see D. 

C) q 2 and g 3 variable, fixed masses 

Instead of holding charges fixed one can fix the masses and study stability in the q 2 , q 3 plane. 
One particular case is m 2 = m 3 = 00 where one has the Born-Oppenheimer limit and one has 
the diagram calculated by Hogreve (Fig. 7). 
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In the g2, qs plane there is again for the general mass case a dividing line where the binding 
energies of the two subsystems (12) and (13) are equal: 

m 2 2 m 3 
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mi + m2 



mi + ^3 



In the two sectors thus defined there are two instability regions for which we have been able to 
derive a new concavity property: 

Theorem IV 

Define z 2 = l/ffe, 23 = l/q^, the image of q 2 > q 3 > is z 2 > z 3 > 0. Then, in the z 
variables the two instability regions are convex (Fig. 8). The proof is based on a rescaling such 
that the binding energy of the relevant subsystem remains constant on a segment in the £2^3 
plane. 



Unstable 


Stable / 




Unstable 



z 2 



Figure 8: 



D) An illustration: the instability of the systems ape or ap/i 

In the Born-Oppenheimer limit it is known that such systems are unstable [H. Spruch and 
collaborators |l5j have given arguments which seem to indicate that this might remain true for 
the actual masses of the protons and of the a particle, but, to our knowledge, a completely 
rigorous proof does not exist. 
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STEP 1 Take m 2 = m 3 . Then if q 2 = q 3 we have instability if g 2 = q 3 > 1.24 and 
m% = m 3 = and m 2 = m 3 = oo, and by concavity for any m 2 = m 3 . Now consider the 
segment q 2 = 1-24, < q 3 < 1.24. Along this segment the subsystem with the most negative 
binding energy is (12) and this energy is constant. There is no stability for q 2 — q 3 = 1.24, and 
no stability for q 2 = 1.24, q 3 very close to zero, because then particle 3 is submitted to a very 
weak force and is therefore very far away most of the time while 12 is overall repulsive for 3. By 
concavity there is no stability for q 2 = 1.24, < g 3 < 1.24. The same kind of argument applies 
to 1.24 < q 3 < oo, because one has instability for q 3 — > oo q 2 > 1, and one can use concavity in 
the inverse charge. The conclusion is that if m 2 = m 3 , one has no stability if either q 2 > 1.24 
or q 3 > 1.24. 

STEP 2 Assume q 2 > q 3 . Then, in the whole sector, m 2 > m 3 , the lowest two-body threshold 
is given by the (12) system, therefore in the a triangle the region a 2 < a 3 , is completely unstable 
for q 2 > 1.24 by use of the star shaped property. 

The systems 

ape~ apfi~ 
ade~ ad\i~ 
ate~ atjjT 

satisfy precisely the conditions: q 2 > 1.24, q 3 — 1, m 2 > m 3 , and are therefore unstable, in 
the sense we have given to "instability" . Notice that the proof would fail if the particle with 
charge 2 was lighter than the particle with charge 1. 

However, as pointed out by for instance Gerstein [ 16| , some of the levels of these systems 



are "quasi stable" in the Born-Oppenheimer approximation in the sense that the minimum of 
the Born-Oppenheimer potential is below the value it takes for infinite separation between the 
two nuclei where one of the limit atomic states is excited and degenerate with the other one. 



Gerstein [17] went as far as estimating the lifetimes of these quasi-stable states and showed 



that the lifetime increases drastically when the proton is replaced by a triton. One should also 



mention the metastability, where the Born-Oppenheimer curve has a minimum above zero [12 



4 Four-body case: equal charges 



Most of what I will say concerns systems of two positive and two negative charges of absolute 
value e. 

However, let me start with the case 

pe~e~e~ , 

i.e., a doubly negative hydrogen ion. 



Such a state according to a review by Hogreve [|18[ does not seem to exist. In the limit of 
an infinitely heavy proton, the Lieb bound on n, the number of electrons around a charge Z 
DJ, n < 2Z + 1, which is a strict inequality, gives n < 3. In fact no doubly negative atomic 
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ions seem to exist in nature, while singly negative ions may (like H ) or may not exist (like 
the case of the rare gases) . 

We return now to systems with charges — e — e + e + e, and first of all m~m~ M + M + , i.e., two 
negatively charged particles with equal mass and two positively charged particles with equal 
masses. A familiar example is the hydrogen molecule e~e~p + p + . A more exotic example is the 
positronium molecule e~e~e + e + . 

It has been realized by Jurg Frohlich that up to very recently there did not exist any 
rigorous proof of the stability of the hydrogen molecule. It was believed to be stable because 
of experiment of course, and of Born-Oppenheimer calculations. 

Two groups (Frohlich et al., Richard) investigated this problem and finally joined their 
efforts to produce a completely rigorous proof fl9 |. 



The simplest approach, whose idea comes from J.-M. Richard [2(| consists of starting from 
the work of 0re, which is valid by scaling for a system A~A~A + A~^ 
variational trial function, of the form 



21[ . 0re used a very simple 



V> = exp -i (r 13 + n 4 + r 23 + r 2 ^j 



cosh 



ri3 - r M - r 23 + r 24 



Notice that the distances between particles with same charge sign do not appear. All 
integrals can be carried analytically and it is found that the energy is less than 



2.0168 E (A + A~ 



The system is therefore stable because it cannot dissociate into A + A~ + A + A + . It cannot 
dissociate either in A + A~A + + A~, because between these two systems there is a long-distance 
Coulomb force, producing unavoidably infinitely many bound states. 

If we take now 

1 1 2 

m e m p rriA 

we see that the binding energy of e~p is the same as that of A + A~. However, 

' I ' I i • I iy iy I ry I ■ I rp 

tAJ p | J^g Jb p | •-t' g p I J^p I wL'g 



E\X e ^ X e , Xp, %p) EyXp^ Xp^ X e , X e J <^ E 

by concavity in the inverse masses. 
So, 

E(p + ,p + ,e~,e-) < E(A + ,A + ,A~,A-) < 2E{A + A~ 
Hence, ppee is stable. 
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One can wonder if stability remains if the masses of two particles of the same charge are 
different, i.e., 

A + B + C-C- 

Then there is still a unique possible dissociation threshold: 

A + C- + B + C- 

0re has predicted explicitly that the system pe + e~~e~ is stable [f22| , and this has been observed 
experimentally by Schrader and collaborators E3|. 



It is also easy, from the upper bound of the energy of e + e + e e to show that the systems 



p d e e , p t e e , d t e e are stable. This is implicit in the work of Richard p0|| , established 
in the thesis of Seifert and I present here my own version. By concavity we have 



E(x A ,x B ,xc,x c ) < E 



XA+XB XA+XB 



,x c ,x c 



fPI XA+XB I XQ_ XA+XB I XQ_ XA+XB I X£_ XA+XB I Xc 

^ 1 4 ~ r 2' 4 ~ l ~2' 4 ~ r 2' 4 ~ r 2 



< -2.0168 | 



If the inequality 



-2.0168 



X'A + X-fl _|_ £^ 

4 ~ r 2 



< - 



XA , X£ 

. 2 2 



+ 



2 2 



is satisfied, the system is stable. If > rug > mc one finds that this condition is certainly 
satisfied if m B > 5mc- 



Using the more refined bound |19 



Ea+a+a-a- < -2.06392E{A + A + ) , 

which uses more sophisticated trial function and must be "cleaned" from numerical roundup 
errors, one gets 

tub > 2.45mc . 



However, Varga and collaborators p5| , using trial functions leading to integrals which can 
be expressed analytically, and adjusting parameters, have found that one has stability for any 
rriA and tub, including the case where one or two of them are less than mc- By a tedious but 
feasible exercize, one could, using concavity, transform this calculation, which is unavoidably 
done for discrete values of the masses, into a very inelegant proof. Let us hope that someone, 
in the future, will find a still more clever trial function and avoid this. 

Let us mention, finally, a conjecture by Jean-Marc Richard: "A four-body system with 
two positive-charge and two negative-charge particles (equal in absolute value) is stable if one 
three-body subsystem is stable" . 
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